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Abstract 

>Y 

■ The Cauchy problem is considered for the scalar wave equation in the Kerr geome- 

try. We prove that by choosing a suitable wave packet as initial data, one can extract 
energy from the black hole, thereby putting supperradiance, the wave analogue of the 
Penrose process, into a rigorous mathematical framework. We quantify the maximal 
energy gain. We also compute the infinitesimal change of mass and angular momentum 
of the black hole, in agreement with Christodoulou's result for the Penrose process. 
The main mathematical tool is our previously derived integral representation of the 
00 ' wave propagator. 
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1 Introduction and Statement of Results 

Near a rotating black hole there is the counter-intuitive effect that the physical energy of 
"q I particles or waves need not be positive. As discovered by Roger Penrose [T2], this effect 

can be used to extract energy from the black hole. In the process proposed by Penrose, a 
classical particle enters the so-called ergosphere, a region outside the event horizon, where 
it disintegrates into two particles. One particle of negative energy falls into the black hole, 
whereas the other particle has energy greater than the initial energy, and escapes to infinity. 
^ The wave analogue of the Penrose process is called superradiance; it was proposed in [16] 

and first studied in [HI [14] and [UH], see also jSJQTJQl)]. In these papers, superradiance 
was considered only on the level of modes, i.e., by considering the reflection coefficients for 
the radial ODE obtained after separating out the time and angular dependence. A more 
convincing treatment of superradiance is to consider the Cauchy problem for "wave packet" 
initial data, and in [TJ, this was carried out numerically for the scalar wave equation. In 
this paper we shall consider the Cauchy problem analytically, giving a mathematically 
rigorous treatment of superradiance for scalar waves. Our analysis is based on the integral 
representation for the wave propagator obtained in [U [9] . 

A rotating black hole is modeled by the Kerr metric, which in Boyer-Lindquist coor- 
dinates (i, r, if) with r > 0, < $ < 7r, < (p < 2ir, takes the form 

ds 2 = gjk dx^ dx k 
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A ,1. ■ o „ , ,9 TT /dr 2 . , q2 \ sin 2 d 2 2 



_ 7 (eft - a sin 2 d di^) 2 - U ^— + dd 2 J — (a (it - (r 2 + a 2 ) dp) . (1.1) 

Here M > and aM > denote the mass and the angular momentum of the black hole, 
respectively, and the functions U and A are given by 

U(r,fi) = r 2 + a 2 cos 2 d, A(r) = r 2 - 2Mr + a 2 . 

We consider only the non-extreme case M 2 > a 2 , where the function A has two distinct 
roots. The largest root 

n = M + VM 2 - a 2 

defines the the event horizon of the black hole. We restrict attention to the region r > T\ 
outside the event horizon where A > 0. The metric (jl.ip does not depend on t nor ip and 
is thus stationary and axisymmetric, admitting the two commuting Killing fields and 
tj^. The ergosphere is defined to be the region where the Killing field is space- like, that 
is where 

r 2 - 2Mr + a 2 cos 2 d < . (1.2) 

The ergosphere lies outside the event horizon r = T\ , and its boundary intersects the event 
horizon at the poles d = 0, ir. 

We now briefly review the Penrose process (for more details see [15] ) . The 4-momentum 
pi of a massive point particle is time-like and future-directed, and thus its energy (p, Jj) is 
positive whenever Jr is time-like. However, in the ergosphere the vector becomes space- 
like, and hence the energy of the point particle need not be positive. Penrose considers 
a particle of positive energy which splits inside the ergosphere into two particles whose 
energies have opposite signs. By finely tuning the energy and momenta of these particles, 
one can arrange that the particle of negative energy crosses the event horizon and reduces 
the angular momentum of the black hole, whereas the other particle escapes to infinity, 
carrying (due to energy conservation) more energy than the original particle. In this 
way, one can extract energy from the black hole, at the cost of reducing its angular 
momentum. Christodoulou [3j showed that the infinitesimal changes of mass 5M and 
angular momentum 5{aM) of the black hole satisfy the inequalities 

2 2 

SUM) < ri+a 5M < 0, (1.3) 
a 

and as a consequence he showed that it is not possible to reduce the mass of the black 
hole via the Penrose process below the irreducible mass M irr defined by 

M 2 r = i (M 2 + ^M 4 - (aM) 2 ) . (1.4) 

We now recall the wave equation in the Kerr geometry and its separation; for more 
details see [8]. The scalar wave equation is 

d ( / _ — i;j d 
V— det g dx % \ dxi 

and in Boyer-Lindquist coordinates this becomes 



^' V 4 V^ = ^4=£-(^^^^)cI> = 0, (1.5) 



d . d 1//9 <a<9 <9\ 2 d , 9 _ d 
A— + — \(r 2 + a 2 )— +a—\ - sm 2 d 



dr dr Ay dt dip J d cos d d cos d 

1 ( ■ 2 q d d 
asm v— — h — — 



sin 2 d V dt dip 



$ = 0. (1.6) 
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Using the ansatz 

$(t, r, <p) = e-^-^ R(r) 0(0) (1.7) 

with u; € K and fc G Z, the wave equation can be separated into both angular and radial 
ODEs, 

K^ k Rx = -\Rx, Au,k @A = A ©a • (1.8) 

The angular operator .A^.fc is also called the spheroidal wave operator. It has a purely 
discrete spectrum of non-degenerate eigenvalues < Ai < A2,.... The corresponding 
eigenfunctions 0«' fc are referred to as spheroidal wave functions. In order to bring the 
radial equation into a convenient form, we introduce a new radial function cj)(r) by 

<p(r) = yV 2 + a 2 R(r) , (1.9) 
and define the Regge- Wheeler variable u S K by 

du r 2 + a 2 



dr A 



(1.10) 



mapping the event horizon to u = —00. The radial equation then takes the form of the 
Schrodinger equation 

<h +V(u)) <j>(u) = (1.11) 



du 2 



with the potential 



= " (" + -^t) + t¥^4 + ~4== % dl^oJ . (1.12) 

V r z + a z J (r z + a z ) z ^ r 2 + a 2 

Starobinsky [13J analyzes superradiance on the level of modes. Using the notation in [9], 
we can explain his ideas and results as follows. We fix the separation constants k > 0, u> 
and A ra . Introducing the notation 

ft = u) — u>o with ujq = ^ k , (1-13) 

rf + a z 

the potential (|1.12p has the asymptotics 

lim V(u) = -n 2 , lim V(u) = -to 2 . 



Thus there are fundamental solutions <f) and 4> of (jl.lip which behave asymptotically like 
plane waves, 

lim e- inu 4>(u) = 1, lim (e~ mu ku))' = (1.14) 

u— >— 00 u— >— 00 V / 

lim e iuJU ${u) = 1, lim (e^ u ku))' = 0, (1.15) 

(see [9j for details). The solution <f> = <j) has, near the event horizon u = —00, the 
asymptotics (p(u) = e - *^". Taking into account the factor e~ lujt in the separation, this 
corresponds to a plane wave entering the black hole. A general solution can be expressed 
as a linear combination of <fi and its complex conjugate, 

(j){u) = A s (j){u)+B^{u) . (1.16) 
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Equation (|1.15p shows that the solution cj) behaves near infinity like <f){u) = e luJu . Hence 
the corresponding time-dependent solution behaves like the plane wave e - luJ ( t + u ) ari d cor- 
responds to an incoming wave propagating from infinity. Likewise, the solution <j) corre- 
sponds to an outgoing wave propagating towards infinity. Computing the Wronskian of (j) 
and 4> near the event horizon and near infinity gives the relation 

\A\ 2 -\B\ 2 = -. (1.17) 
u 

The quantities u; 2 |A| 2 and w 2 |i?| 2 have the interpretations as the energy flux of the incom- 
ing and outgoing waves, respectively. If the right side of (|1.17p is positive, the outgoing flux 
is smaller than the incoming flux, and this corresponds to ordinary scattering. However, if 
the right side of (|1.17p is negative, then the outgoing flux is larger than the incoming flux; 
this is termed superradiant scattering. According to (|1.17p . superradiant scattering ap- 
pears precisely when u and f2 have opposite signs. Using (|1.13p . we see that superradiant 
scattering occurs if and only if u lies in one of the following intervals 

ujq < uj < , < u < ujq , (1-18) 

depending on the sign of ujq. The gain in energy is determined by the quotient of outgoing 
and incoming flux, 

* = w ■ (L19) 

Starobinsky [13] made an asymptotic expansion for 91 and found a relative gain of energy 
of about 5% for k = 1 and less than 1% for k > 2. This was verified later numerically 
in pp. Teukolsky and Press [H] made a similar mode analysis for higher spin and found 
numerically an energy gain of at most 4.4% for Maxwell (s = 1) and up to 138% for 
gravitational waves (s = 2). 

Our main result makes the above mode argument for the scalar wave equation mathe- 
matically precise in the setting of the Cauchy problem. To state our result, we combine $ 
and dt<&, as in [8], into a two-component vector ^ = (<3?,z9j$). We always assume that 
the initial data is smooth and compactly supported outside the event horizon, 

*o = (*,i$*)|t=o G Wi,»)x5 2 ) 2 . 

The physical energy of the wave is then given by <^, ^>, where <., .> is the energy scalar 
product 

roo rl ( / (~2 I „2\2 \ 

= J dr J d(cos 0) | i^- — - a 2 sin 2 #J d t <S> dt& 

+A-i\$d r & + sin 2 ^^Scos^' + (-r^ ~ ^) dp$d ip &\ . (1.20) 
Provided that the limit t — > oo exists, the energy radiated to infinity can be defined by 

E out = lim <tf(t),X(2n,oo)(r) *(*)>. (1-21) 

where x is the characteristic function. Note that, according to the pointwise decay result 
in [9], we could replace the argument 2r\ by any radius greater than r%. Moreover, again 
using pointwise local decay, the boundary term at r = 2r\ which arises when differentiating 
the characteristic function X(2n,oo); vanishes in the limit t — * oo. We can now state our 
main theorem. 
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Theorem 1.1 Consider the Cauchy problem for the scalar wave equation in the non- 
extreme Kerr geometry for initial data with compact support outside the event horizon, 



* = (*,i$*)ii=o e c^({ ri ,^)xs 2 y , 

having energy <^q,^q>. Fix k > 0, n £ N and u £ (loq,0). Then for any R > r\ 

and 5 > there is initial data \&o G Cq°((R, oo) x S 2 ) 2 suc/i i/iai i/ie Zimit in rti.ffi)) exists 
and 

Eout 



< 6 



<*o,*o> 

with y{ as in U.19\) . The same inequality holds for k < and a; £ (0,luq). 

We emphasize that we allow the initial data to be supported arbitrarily far away from the 
event horizon. This is important in order to avoid artificial initial data which would not 
correspond to an energy extraction mechanism. For example, if one allows the support 
of the initial data to intersect the ergosphere, one could take initial data close to a wave 
packet with zero total energy, in which case the quotient E out /<^ , ^ > could be made 
arbitrarily large. 

We also compute the energy E^ and the angular momentum of the wave compo- 
nent crossing the event horizon of the black hole to obtain the following theorem. 

Theorem 1.2 For any initial data G Cg°((ri,oo) x S 2 ) 2 , the quantities E^ and 
defined by \2.8l \ 5.1\) satisfy the inequality 

r 2 -\- a 2 
a 



This shows explicitly that Christodoulou's estimates (|1.3l II. 4ft also hold for energy extrac- 
tion by scalar waves, in agreement with Hawking's area theorem |15j . 



2 Absorbtion of Energy by the Black Hole 

Recall from [9 J that given initial data G C^°(M x S 2 ) 2 , the solution of the Cauchy 
problem for the scalar wave equation (jl.6p can be represented as 

k£l neTN J ~°° a,b=l 

where the sums and the integrals converge in L 2 oc . In the above integral representation, 
the coefficients t a b are defined by 

OL OL OL 

t n = 1 + Re-, t 12 = t 21 = -Im-, *22 = 1-Re-, (2.2) 
and the complex coefficients a and (3 are defined by 

= a<j> + P<j>, (2.3) 

where <j>(u) and 4>(u) are the fundamental solutions of the radial equation having the 
asymptotics (|1.14p and (|1.15p . respectively. We will refer to the complex coefficients a 
and (3 as transmission coefficients. Finally the functions ^? wn (r, a = 1,2 are the 
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solutions of the wave equation (|1.6p . with fixed quantum numbers k,u),n, corresponding 
to the real-valued fundamental solutions of the radial equation given by 

(p 1 = Re(p , (p 2 =Im(p. (2.4) 

Here we shall always consider a fixed fc-mode, and without loss of generality we assume 
that k > 0. For notational simplicity, from now on we omit the index k and the in- 
dependence. Furthermore, we introduce the vector notation 

*> - ( f^; Jj ) , vr == ( tH ) . t - (t- w . m 

which allows us to write the integral representation (|2.1|) in the compact form 

*(t, r, 0) = — / "o e"^(*" n (r, 0), T ) ca . (2.6) 



We now introduce some basic quantities needed for the formulation of the superradi- 

"fOC 





ance property for wave packets. The total energy E tot of an initial data £ Cq°(R x S" 2 ) 2 



for the scalar wave equation (jl.6p is defined by 

^tot = <*o,*o> • (2.7) 

This energy is conserved throughout the time evolution of the scalar wave, meaning that 
if ^(t) denotes the solution of the Cauchy problem for the scalar wave equation, then 

<* 0! #o> = <#(i),#(t)> • 

Recall that the outgoing energy, which represents the energy radiated to infinity, has been 
defined in (|1.2ip . Finally, the energy absorbed by the black hole is defined by 



E 



bh 



lim <tt(t),x ( _ 0Q ,ar 1 ) ( r )*(*)> 

t— >oo v ' ' 



We next derive useful expressions for these quantities in terms of the initial data and 
the transmission coefficients a and (3. The expression for 2£ to t is an immediate consequence 

of dZSD - 

Proposition 2.1 Choose a fixed k G Z suc/i f/iai u>o < 0. T/ien 

n.GlN J 00 

where the series converges absolutely. 

We next want to compute E^h and 2£ ut- First, we need an argument which ensures 
that the sum over the angular momentum modes converges. Since E out = E to t — -^bh 
and the convergence of the angular sum is already known for the total energy according 
to Proposition 12.14 it suffices to consider for example E ont . In [lOj it is shown that the 
outgoing energy is bounded uniformly in time, i.e. 

<*(t),X(2n,oo)(r)*(*)> < far all t. 
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Moreover, it is also shown in [TU] that the outgoing wave can be approximated by a finite 
number of angular momentum modes uniformly in time, i.e. for any 5 > there is N such 
that 

<* iV (t),X(2r 1 ,oo)(r)* JV (i)> < 6 foralH, 

where 

1 roo J 2 

9 N (t, r,0,<p) = — Y, / 7^ E C n »T (r, *) <«T, *o> • 

n>N a,b=l 

This estimate will always allow us in what follows to interchange the limit t — > oo with 
the sum over the angular momentum modes. 

To compute 2£bh> we make use of the following lemma, which generalizes [3 Lemma 9.1]. 



Lemma 2.2 For any Schwartz function f £ S(R x R) and any uq £ M., we define the four 
functions A s ^ s / , where s,s' £ { — 1, 1}, by 

i*uq /'OO /'OO 

A s>s , = lim / du I duj / duj'e-^-^ 1 er isnu+is ' n ' u f(uj,uj') . 



OO J — OO 



Then 

A s , s , = 2n(6l6A f(u,u)tb;, 

^ ' J —OO 

where 5f denotes the Kronecker delta. 

Proof. Noting that VL — Q' = uj — u/, the case s = s' was proved in Lemma 9.1]. Thus 
it remains to show that 

rUQ roo roo 

lim / du I du I du'e-^-"'* e ±i( - n+n >f(uj,uj') = 0. 



OO J — OO J —OO 



Noting 

e -i(u-w')t e ±i(n+n')u _ 1 rQ _|_ i\(g , _|_ ^\ e -i{w-ui')t±i(n+n')u 

(— it ± iu + 1) (it ± m + 1) w w 

we can integrate the u>, u/-derivatives by parts. Applying Fubini's theorem and estimating 
the resulting ^-integral for t > 2\uq\, we find 

7— ,4. TC *< n+n >du = / ((t + u) 2 + lp((i-«) 2 + lp d« 

{tTu + i){t±u-i) 

1 f" -1 

< Vet'* ((t + uf + i) 4 du. 

J-oo 

Since the last integral is bounded uniformly in t, the entire expression tends to zero as 
t — ► oo. ■ 
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Proposition 2.3 Choose a fixed integer k > 0. Then the limit in 112. ty) converges and 

1 fOO J 



;LUn\ 

/C 2 



n.eIN 

where P is £/ie matrix 

2 V -* 1 

T/ie above series is absolutely convergent. 

Proof. Substituting the integral representation (12. ip into (12. 8ft . we find 

1 /' 2ri , /*°° du' f°° dtu j(u] n, - u , . 



where 



n,n'6Na,6,c,d=l J s2 



b i / J 



and £{^f% n , denotes the bilinear form corresponding to the energy density as given 
in (H3D]). 

We next justify that we may interchange the limit t — * oo with the sums over the 
angular momentum modes and that the series converge absolutely. According to [10} 
Theorem 8.1], the energy of the large angular momentum modes near infinity is small 
uniformly in time. Due to conservation of energy and the local decay [9], a similar result 
follows for the energy near the event horizon. Hence the integrals in (|2,9p can be approx- 
imated uniformly in time by a finite number of angular momentum modes. Thus in what 
follows we may restrict attention to fixed angular momentum modes n and n' . 

Near the event horizon, we can expand the coefficient functions in the energy den- 
sity £(*£ v , to obtain 



r z _|_ a z 

a 2 k 2 



= (r? + a 2 ) *™ + d u ^' n d u ^ n ) ~ -f^2 C + 0^'(e 7 ") 

= 8„ v e m ^^pr+wp^r - (r ^ fc Q 2 2)2 f r) + ew(e 7n ) 

with the fixed constant 7 > as in [9l Eq. (3.9)]. Here the error term can be bounded by 



CW(e 7 ")l < Ce lu (\r b ' n '\ + W4>t' n '\ + \d u <Pt' n '\) {\<PT\ + Mrl + \du<t>T\) (2-10) 



with a constant C independent of uj and a/. Using the Jost estimates of [9, Theorem 3.1], 
the fundamental solutions have the following expansion near the event horizon, 

= c os(nu) + O w y(e^ u ), = sin(ftu) + 6^,{e~< u ) . 
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Here the error term is bounded by 

\0^(e^ u )\ < Ce^, \d u 6^,(e^)\ < C (1 + \oj\)e? u , 

where the constant C is again uniform in to (see (9J eqns (3.7, 3.10)]). This estimate also 
shows that the error term (|2.10p is bounded by 

|CW(e 7n )| < Ce^ (1 + |a/|)(l + M) . (2.11) 

We thus conclude that 
A 



r 2 + a 2 



where O UJ ^i{e lu ) satisfies (|2.11[) (possibly with a new constant C) and 

ipi(u) = cos(Slu) , ip2 ( u ) = sm(Qu) . 

Using the asymptotic formula f)2 . 12|) in (|2,9p and writing the radial integral in terms of the 
Regge- Wheeler variable u, the factor A/(r 2 + a 2 ) drops out. Applying [lQl Lemma 3.1], 
we can write E^h as 

fUQ poo roc 

E hh = lim / du / duj' / dwe - ^^ 



^bh 

t— >oo 



2 



with Schwartz functions /b iC and 5. To treat the error term, we first apply Fubini's theorem 
to exchange the orders of integration. This gives 



lim / dJ / du e-^-^ g(u,u') 



and using (|2.1ip . we can apply the Riemann-Lebesgue to conclude that this limit is zero. 
The term involving the f^ c can be treated by applying Lemma [2.21 Using the orthogonality 
of the spheroidal wave functions, we obtain 

1 f°° dco ^ /TM1 TPT TUn . f 2 , a 2 k 2 



(r? + a 2 ) 2 



We finally apply the identity 



n 2 k 2 

2,r>2 a 2,r>2 2 2,2 on 

ijJ + S Z — —5 — — = UJ + \ l — OJn = UJ + CO — OJOJq = ZUi I . 

(rf + a 2 ) 2 



Using the matrix identity 

T - TPT = TQT with Q := T" 1 - P , 
a short calculation yields the following result. 
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Proposition 2.4 Choose a fixed fcGZ such that uq < 0. Then the limit t — > oo m 111.21]) 

exists and 



neiN 

where Q is i/ie matrix 



1 />oo J 

^ ^-TKT ■/ — OO 



f |a-/3| 2 i(a + /3) (a - /?) 



2w V -*(«- P) ( a + 0) |a + /?| 2 
3 Energy Propagation of Wave Packets near Infinity 

We fix u> £ (u}q,0) and h E N and consider initial data V&o i n the form of a linear combi- 
nation of outgoing and incoming wave packets, 



*o = 
where L > 



id t <5> 



- ft. .(, 9 ) ™ . 



! ) +cbute* 5u ( *~ 
a; / V — <jj 



(3.1) 



1 fu-L 2 

T)l(U) = —j=7] 



L \ L 



with r/ G Cq°((— 1, 1)) being a smooth cut-off function, and where 6^5(1?) is an eigenfunc- 
tion of the angular operator A. 

In the next proposition we compute E tot asymptotically as L —* 00. 

Proposition 3.1 For £/ie initial data given by the wave packet \3.1\) . 

lim E tot (*o) = -7- £ f lim / — (*o n >Dc 2 , (3-2) 

L— >oo 47T ^ — * L^oo l-~U 

neIN 

*r - f^'^V (3.3) 



<^„,^o> 



TTie above series converges absolutely. 



Proof. Applying |10| Lemma 8.3] with / = 1, we see that the energy of the large angular 
momentum modes is small uniformly in L. This allows us to interchange the limit L — » 00 
with the sum over the angular momentum modes. Hence in what follows we may restrict 
attention to a fixed angular momentum mode n. 

Since the wave packet is localized near infinity, it is preferable to work with the funda- 
mental solutions <fi having the asymptotics (|1.15j) . They are related to the functions (fr-m 
(gaD by 

(iMS) 

where 

a + (3 i(a — (3) 



A :-- 



a + j3 —i(a — (3) 
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In order to compute the inverse of A, we first note that, from the asymptotics (|1.14|) 
and (|1.15p of the fundamental solutions of the radial equation at infinity and at the event 
horizon, we obtain the following Wronskians, 

w(4>,4>) = —2iu>, w((p,(j)) = 2iQ, . 
On the other hand, from (|2.3p . we know that 

= (\a\ 2 -\p\ 2 )w(lJ) . 

The last two identities imply that 

H 2 -|/3| 2 = (3-5) 

Using this relation, we obtain 



B : = A" 1 = — ^ ~® + f3 ~ a + P 



2uj \ i(a + (3) -i(a + (3) 
Using again the notation (|2.5p . we have 

= B*£ n , (3.6) 

so we can write the total energy as 

i x-^ r°° dto v f 

ng]N J 00 

The quadratic form can be written in components as 



'o /C 2 



(^B'TB = (*^) 1 (B t TB) 11 (*^) 1 + (*-«) 2 (B t TB) 22 (*D 2 (3.7) 

+(i^ r )i(B t TB)i 2 (*r)2 + (i|r)2(B*TB) 21 (*ni . (3.8) 

and the matrix B*TB is given by 

_ to f 1 -a/p 



BTB , , 

2uj V -<x/P 1 



We consider the contribution by the diagonal terms (|3.7j) and the off-diagonal terms (J3T 
separately. Re-expressing the diagonal elements in terms of ^q 11 we obtain 



and the corresponding contribution to E tot is given by (|3.2p . 

It remains to show that the contribution of the off-diagonal terms (|3.8|) to E'tot tends 
to zero as L — > oo. Intuitively, this can be understood from the fact that the off-diagonal 
terms involve oscillatory terms e zizlu}L , whose integral becomes small for large L. In 
order to make the argument precise, we use an argument based on the Riemann-Lebesgue 
Lemma. We only consider the (1, 2) element; the argument for the (2, 1) element is similar. 
We write 



1 



-- (*-") 1 (B t TB) 12 (*r)2 = e- 2i " L2 g L {u) 
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with 



0iM = ~j <^o,e lujL2 ^ n >^j (± <e^ 2 ^«,* >) • (3.9) 

For clarity we point out that the L-dependence of g^ comes from \&o (see (|3.ip ). as well 
as from the phase factor e tLuL . The corresponding contribution to E tot is obtained by 
integrating over u, and thus the remaining task is to show that 

/oo 
e~ 2iujL2 g L {u) dw = 0. 
-oo 

This is a direct consequence of the following two lemmas. ■ 



Lemma 3.2 Let Jl € -L 1 (R), L > 0, be a family of functions with which converges in 
L^R) as L — > oo. Then 

/oo 
e~ 2iujL2 f L (u) du = 0. 
-oo 

Proof. Let / = limi,_ +00 in L . Using the inequality 

/oo 
\f - f L \(u) du , 
-oo 

the first term tends to zero by the Riemann-Lebesgue lemma, while the second term tends 
to zero by hypothesis. ■ 



e" 2 ^ f L (u) du 



< 



-2iu)L 



f{uj) du 



Lemma 3.3 For any n,n and Co, the functions g^ converge in L 1 (M) as L — > cxd. 

Proof. In view of Lebesgue's dominated convergence theorem, it suffices to show that 

(a) The pointwise limit limx-^oo ffL^) exists for all co ^ 0, and 

(b) There is a constant C such that for all u and L, 

\9M\ < . 

To show (a), note that in the limit L — > oo, the support of moves to infinity. Thus using 
the plane wave asymptotics of the Jost solutions ip^™ together with the error estimates 
in the proof of Lemma 3.3 in [9], a straightforward calculation shows that both brackets 
in (|3.9p converge for any fixed to ^ 0. This proves (a). 

To prove (b), we first analyze the behavior near oj = 0. First note that the factors lo^ 1 
in the brackets in f|3.9|) are compensated because the energy scalar product involves a 
factor of u (see [HI Eq. (2.14)]). Thus the estimates for the Jost functions in Lemma 3.6 
in [9] yield that the two brackets in (|3.9[) are bounded for small |o>|, uniformly in L. Next, 
the convexity argument in [9j Section 5] yields that the factor a/ (3 is bounded near u = 0. 

On any compact set which does not contain u = 0, we know from [9l Section 7] that the 
factor a//? is a continuous function. Since the Jost solutions form holomorphic families and 
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the estimates of Lemma 3.3] hold, it follows that the integrand in (|3.9p is continuous 
and thus bounded on any compact set, uniformly in L. 

To control the behavior for large \oj\, we first see from (|3 . 5[) that \a//3\ < 1. Thus 
to finish the proof of (b) it suffices to show that the energy scalar products in (|3.9p are 
bounded, uniformly in L. Indeed, we shall show that they have rapid decay in u, uniformly 
in L, i.e. for any p£N there is a constant c p such that 



<* ,^ n > 



<f un , ^ > 



< 



\uj\P 



VL, Vw with |u;| > 1 . 



(3.10) 



For this we modify the argument in the proof of |1Q|. Lemma 3.1 (iii)]. We write the radial 

ode (frrni as 

Iterating this relation gives 

u; 21 ^ = (-d 2 u + {V + u 2 )) 1 ^. 
We thus obtain for suitable functions F and G, which are independent of cj, the formula 

w 2i $™ = (-flg + cjF + G) 1 ^>" n , 
: i$$<"*) and is defined by $L7\ ES}. Hence 



where ^ u 



1 



LO 



LO 



21 



<{(-dl + ujF + Gyyy ,^ n > , 



(3.11) 



where star denotes the formal adjoint obtained by partial integration. Writing the func- 
tion ((— d 2 +LoF+G)*) l ^Q as a polynomial in u, each coefficient is again a smooth function 
in u, which is bounded by const/\/L and is supported in the interval u G [L 2 — L, L 2 + L\. 
Hence again using the Jost estimates in [9j Section 3.1], we conclude that 



<((-flg + ojF + G)*) l * , * wn > 



where the constant C depends only on \&o and I, but is independent of L. Since I is 
arbitrary, we obtain the desired estimate for the first summand in (13,101) , The proof for 
the second argument is similar. ■ 



For later use, we restate the inequality (|3.10p in terms of $™ as defined by (13. 3D . 

Remark 3.4 For any n, h E N and Co G (loo,0), the function $™ has rapid decay in lu, 
uniformly in L, i.e. for any p G N there is a constant c p such that for all u and all L, 



< 



The result of Proposition 13.11 can be understood more directly from the following 
intuitive considerations. For large L the wave packet is localized near infinity, and thus 
its energy is well-approximated by the energy in Minkowski space. Considering only one 
angular mode and integrating out the angular variables, we thus obtain 



tot 



\dtM 2 + \drM 2 ) du + OiL- 1 ), 
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(where the angular derivatives are contained in the error term). Writing the solution of 
the one-dimensional wave equation as a Fourier integral involving left- and right-going 
waves, 

<Kt, u) = ^ e** (J>l(u;) e~^ u + $ R (u) e*"*) du , 
we can rewrite E+~t. as 



tot 

oo 



E tot = i y°° u? (\h(u)\ 2 + |^H| 2 ) dw + OiL- 1 ) . (3.12) 

Next, using the asymptotics (|1.15p of the fundamental solutions 4> together with the asymp- 
totic form of the energy density in (|1.20j) . we find that 

h{") = ^2 <* wn .*o> + 0(L- 1 ) , fl ( w ) = -L <¥™ * > + 0(L- 1 ) . 

Using this formula in (|3.12p . we again get the expression in Proposition 13.11 

Using the same method as in Proposition 13. 1\ we next compute the outgoing en- 
ergy E out for our wave packet. 

Proposition 3.5 For the initial data given by the wave packet \3.1\) . 

i pea j 

lim £ out (*o) = -r- £ lim / — (*r,R*o") C2 , (3-13) 

L— >oo 47T *■ — * L^oo l~,LO* 

nglN J ~°° 

where R is the matrix 



lal 2 







R = | I . (3.14) 

1 

The above series converges absolutely. 

Proof. Applying [10} Theorem 8.1], we see that the energy near infinity of the large 
angular momentum modes is small uniformly in L. This allows us to interchange the 
limit L — > oo with the sum over the angular momentum modes. Hence in what follows we 
may restrict attention to a fixed angular momentum mode n. 
Using (|3.6p . we can write the result of Proposition 12.41 as 



I >>oo J 

£°ut = «- £ / -o <*o n .B*TQTB*£ 



ujn\ 

C 2 



ng]N 

The matrix B*TQTB can be computed to be 

n (\a\ 2 /\P\ 2 -a/P 



B TQTB l - 

2a; \ a/p 1 

The diagonal terms give rise to f)3. 13j) . The off-diagonal terms, on the other hand, are 
exactly the same as those considered in Proposition 13.11 Thus we already know that they 
vanish in the limit L — > oo. ■ 
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4 Computation of Energy Gain 



Combining Propositions 13.11 and 13.51 we obtain 

/oo 
-oo 



, . E out 
iim — — 

L^oo E tot 



alto 
C 2 — 2 



Elim 



^2 



(4.1) 



To discuss this formula, let us see how to recover Starobinsky's result (II. 19ft . To this end, 
we replace the superposition in the numerator and denominator by a single wave mode 
with quantum numbers Co ^ and n. Then (|4.1j) simplifies to 



£ out _ (*r>R(^^)*o n )c 2 



tot 



/C 2 



If we consider an outgoing wave at infinity, the vector $™ vanishes in the first component. 
According to (|3.14p . the above quotient gives one. This result is immediately clear, because 
if we take an outgoing wave near infinity, it will not interact with the black hole and simply 
escape to infinity. It is more interesting to consider an incoming wave at infinity. In this 
case, the second component of the vector f ™ vanishes, and thus 



E, 



out 



E, 



tot 



a(Co) 



0{u) 



Expressing A and B in terms of our transmission coefficients a and (3 (by a straightfor- 
ward calculation using (jl.161 HOj) together with (|3.5|) ). this is in complete agreement with 
Starobinsky's result [13]. 

The remaining task is to prove that in the limit L — > oo, our wave packets become 
"more and more localized" near u = u and n = h. 

Proof of Theorem We consider the wave packet initial data (13. ip with Co and n 

equal to oo and n in the statement of the theorem. We choose c, 



out 



and Cj, 



From the asymptotic form of the energy density near infinity, it is obvious that the total 
energy of the wave packet (|3.ip has a non-zero limit as L — > oo. Hence the denominator 
in (14. ip is non-zero. Moreover, since both series in (|4.ip converge absolutely, they may be 
approximated by finite sums. Using the identity 



E, 



out 



i 



E n 



<^o,^o> 



and considering a finite number of modes, we see that our task is to show that 

2\ 



n<no 



lim / (*o"> R(w,n) 



*o )c 2 w2 



^7 = 0. 



(4.2) 



It clearly suffices to prove this for one mode. 

In the case n = n, for given e > we choose an open neighborhood I £ of Co such that 



1 



a(Co) 

W) 



< £ 



for all lo G I F 
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Using the asymptotic form of the energy density and of the fundamental solutions, a short 
computation shows that 

lim (*o ft )2 = pointwise in R 
lim = pointwise in R \ I £ 

L— >oo 

(where ( 1 i r g n )2 denotes the second component of ^q" 1 ), and that the functions (>I/g")2 
and \j/g n are dominated in R resp. in R \ I e by Schwartz functions, uniformly in L (see 
Remark l3.4p . Hence applying Lebesgue's dominated convergence theorem, we obtain 



lim / (*t,[R(LO,n) 
L ^°° Ji F \ 



a(£j) 2 



lim / (*f,*f) c2 ^ 



*0 n )c 2 —o < e ^tot 



Since e is arbitrary, this proves (|4.2p for the summand n = h. 

If re 7^ re, the orthogonality of the spheroidal wave functions for cj = u> together with 
the continuity in u> allows us to choose a neighborhood I £ of uj such that 

<0 n^ )0 n^ >52 < £ forallwe/ e . 
Now we can use the same argument as in the case re = re to conclude the proof of f)4.2|) . ■ 



5 Absorbtion of Angular Momentum by the Black Hole 

We first derive the expression for the angular momentum of the scalar wave. We recall 
from [8J that the Lagrangian is given by 

C = -A|d r $| 2 + ^\((r 2 + a 2 )d t + ad^\ 2 

-sin 2 -d\d cos $(p\ 2 I (a sin 2 #<? t + ci,)$| 2 . 

sin v 

This Lagrangian is axisymmetric. Applying Noether's theorem gives rise to the following 
conserved quantity, 

= dr j d{<xx-ff) JJ ^A, 

where A is given by 

The quantity A can be interpreted as the angular momentum of the wave and A as the 
angular momentum density. 
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Similar to (|2.8|) . the angular momentum absorbed by the black hole is defined by 

A hh = lim / dr / d(costf) / — A[$\ , (5.1) 

provided that the limit exists. In the next proposition we compute A^, again for a fixed 
fc-mode. 

Proposition 5.1 Choose a fixed k 6 Z suc/i £/iai wo < 0. T/ten i/ie limii in A5.1\) exists 
and 



/C 2 ! 



neIN 

where the sum converges absolutely. 

Proof. As in the proof of Proposition l3.5l one sees that we may interchange the limit t — > oo 
with the sum over the angular momentum modes, and that the resulting series converges 
absolutely. Thus it remains to consider a fixed angular momentum mode n. Substituting 
the integral representation (12. 1|) into the definition of A^, we obtain 

A hh = lim dr — — e~<^ >* A u (r > 



where 



2 /■ 



n,n'GN a,b,c,d=l 



and A(^>b' n ' , denotes the bilinear form corresponding to the angular momentum den- 
sity (similar to the bilinear form £ (ty? n , \I / £ m ) appearing in the proof of Proposition 12.31) . 
Near the event horizon, we can expand A(^ n , to obtain 



A{^t' n \ *c U ) = 1 ( r ? + ° 2 ) fe (° + n ') *&'"' + 0{e^ 



2 i 2 V ) c y r< 

r z + cr A 

1 



§ e^e™ k(n + of) <^ n ' € n + o(e' u ) 



Now we can proceed exactly as in the proof of Proposition 



Proof of Theorem Without loss of generality we again restrict attention to the 

case k > 0, so that ujq < 0. We set 

The eigenvalues of the matrix TPT are computed to be zero and 1 + |a| 2 /|/3| 2 . Thus p is 
non-negative. It follows that 
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and using (|1.13|) completes the proof. 
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